Quantum Order: a Quantum Entanglement of Many Particles 
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It is pointed out that quantum states, in general, contain a new kind of orders that cannot be 
characterized by symmetry. A concept of quantum order is introduced to describe such orders. 
As two concrete examples, we discussed quantum orders in free fermion systems and in 2D spin- 
1/2 systems. We generalize the Landau's theory for the classical orders and introduce projective 
symmetry group and Fermi surface topology to characterize different quantum orders. 
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Introduction: Symmetry breaking and the associated 
order parameter have been playing a key role in our un- 
derstanding of phases and phase transitions. Jlj, || How- 
ever, recent study of quantum Hall (QH) liquids reveal 
that QH liquids contain a new kind of order - topological 
order - which can not be characterized by broken sym- 
metry and the associated order parameter.]^, |i] Thus 
Landau's symmetry breaking theory for phase and phase 
transition does not apply to QH liquids and a new theory 
was developed to describe the topological orders in QH 
liquids. H 

The reason that the Landau's symmetry breaking the- 
ory does not apply to QH liquids is because the Lan- 
dau's theory was developed for classical statistical sys- 
tems which are described by positive probability distri- 
bution functions of infinite variables. The QH liquids 
are described by their ground state wave functions which 
are complex functions of infinite variables. Thus it is not 
surprising that QH liquids contain addition structures 
(or a new kind of orders) that cannot be described by 
broken symmetries and the Landau's theory. From this 
point of view, we conclude that any quantum states may 
contain new kind of orders that are beyond broken sym- 
metry characterization. Such kind of orders will be called 
quantum order. [|| 

To visualize the distinction between the classical order 
and the quantum order, we may view the classical world 
described by positive probabilities as a "black and white" 
world, while the quantum world described by complex 
wave functions as a "colorful" world. The Landau's the- 
ory based on symmetry principle and order parameters is 
color blind which can only describe classical orders. We 
need to use new theories, such as the theory of topolog- 
ical/quantum orders, to characterize the rich "color" in 
quantum world. We can also view quantum order as a 
description of the pattern of the quantum entanglement 
in a many-body ground state. A special collective excita- 
tion above a quantum ordered state - gauge fluctuations 
- can be viewed as the fluctuations of quantum entangle- 
ment. In contract, the classical order in a crystal just 
describes a static positional pattern, which has no non- 



trivial quantum entanglement. 

Quantum phase transitions and quantum or- 
ders: Classical orders can be studied through classical 
phase transitions. Classical phase transitions are marked 
by singularities in the free energy density /. The free 
energy density can be calculated through the partition 
function: 
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where h{<f>) is the energy density of the classical system 
and V S pace is the volume of space. 

Similarly, to study quantum orders, we need to study 
quantum phase transition at zero temperature T = 0. 
Here the energy density of the ground state play the role 
of free energy density. A singularity in the ground state 
energy density marks a quantum transition. The simi- 
larity between the ground state energy density and the 
free energy density can be seen clearly in the following 
expression of the energy density of the ground state: 
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where C((f>) is the Lagrangian density of the quantum 
system and V space ti me is the volume of space-time. We 
also note that the free energy density becomes the ground 
state energy density at T = 0. Comparing Eq. ([j]) and 
Eq. (||), we see that a classical system is described by a 
path integral of a positive functional, while a quantum 
system is described by a path integral of a complex func- 
tional. This is the real reason why the classical and quan- 
tum orders are different. According to the point of view 
of quantum order, a quantum phase transition, marked 
by a singularity of the path integral of a complex func- 
tional, in general, cannot be characterized by a change 
of symmetry and the associated order parameter. Thus, 
in general, we cannot use the broken symmetry and the 
Ginzburg-Landau theory to describe a continuous quan- 
tum transition. 

Although the above discussion is limited to zero tem- 
perature, the path integrals of some quantum systems 
can be complex even at finite temperatures. Thus the 
above result also apply to quantum systems at finite tem- 
peratures. It is possible that a continuous phase tran- 
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FIG. 1: Two sets of oriented Fermi surfaces in (a) and (b) 
represent two different quantum orders. The two possible 
transition points between the two quantum order (a) and (b) 
are described by the Fermi surfaces (c) and (d) . 



sition of a quantum system also does not involve any 
change of symmetry even at finite temperatures. 

Quantum orders and quantum transitions in 
free fermion systems: Let us consider free fermion 
system with only the translation symmetry and the U(l) 
symmetry from the fermion number conservation. The 
Hamiltonian has a form 



H =Y1 (4*ii c J + h - c ) 



(3) 



(ij) 



with t*- = tji. The ground state is obtained by filling 
every negative energy state with one fermion. In general, 
the system contains several pieces of Fermi surfaces. 

To understand the quantum order in the free fermion 
ground state, we note that the topology of the Fermi sur- 
faces can changes in two ways as we continuously chang- 
ing tij : (a) a Fermi surface shrinks to zero (Fig. |l]d) and 
(b) two Fermi surfaces join (Fig. |l|c). When a Fermi 
surface is about to disappear in a Z?-dimensional system, 
the ground state energy density has a form 
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where the ... represents non-singular contribution and 
the symmetric matrix M is positive (or negative) defi- 
nite. We find that the ground state energy density has 
a singularity at (i = 0: pe = cf/ 2+D ^ 2 Q(p) + where 
8(a; > 0) = 1, 0(x < 0) = 0. When two Fermi surfaces 
are about to join, the singularity is still determined by 
the above equation, but now M has both negative and 
positive eigenvalues. The ground state energy density 
has a singularity pe = cp( 2+D ^ 2 Q(p) + ... when D is 
odd and pe — c^ 2+D ^ 2 log \p\ + ... when D is even. 

We find that the ground state energy density has a sin- 
gularity at fj, = which is exactly the same place where 
the topology of the Fermi surfaces has a change. Thus 
the topology of the oriented Fermi surface is a "quan- 
tum number" that characterizes the quantum order in 
a free fermion system, (see Fig. |]). A change in the 
topology signals a continuous quantum phase transition. 
Lifshitz|Q] also studied critical properties of such quan- 
tum transitions for D — 3. 

Quantum order in spin liquids: Quantum order 
simply represents the quantum entanglement in ground 
state. The ground state wave function of a free fermion 



system has a simple form of Slater determinant. Its quan- 
tum order can be represented by the topology of the 
Fermi surfaces. In this section, we are going to discuss 
quantum orders in spin liquids. The ground state wave 
functions for spin liquids are more complicated and can- 
not be written as a Slater determinant. Thus we need to 
find a new way to characterize the quantum order in spin 
liquids. 

The problem here is that the many-body wave function 
is too complicated to write down and it is hard to study 
them without writing them down. To gain some insights 
into the quantum entanglements in spin liquid states, in 
the following, we will discuss a concrete representation 
of the spin wave function. We will consider a spin-1/2 
system on a 2D square lattice. For such a system, the spin 
wave function can be written as a bosonic wave function 
$>(xx, xn up ), where X{ is the coordinate of the i th up- 
spin and N up is the total number of the up-spin. 

Within the SU(2) slave-boson approach, ||, [| instead 
of writing down the bosonic wave function <5>({iEi}) di- 
rectly, we regard the boson as a bound state of two 
fermions ip\ and ip2 (which will be called spinons) and 
write $ as 



(4) 



where "Fi^ are the wave functions of ipi,2 and have a form 
of Slater determinant. pX , 11 Actually, more general spin 
wave function can be constructed by introducing a mean- 
field Hamiltonian 



>]4(i)rVi 



(5) 



where ip T = (ipi, i/^), t 1 ' 2,3 are the Pauli matrices, and 
Uij are 2x2 complex matrices. The collection (uij, a l T l ) 
is called a mean-field ansatz. For each mean-field ansatz, 
we can obtain a mean-field ground state by filling the 
lowest 2N up energy levels of H mean with the spinons ip: 



| fl v(«i 3 -,O0' r ')\ 

I * mean j 



*meon° T 5 (Vi, Zi, ■■), where y l are the 



coordinates for ipi and Zi for ^2 ■ The physical spin wave 
function $({a:.;}) can now be obtained by performing a 
projection |lj y i = z i = x i : 



\xi) = \ y - z . =Xi ^[ 



> (y i ,z l ) (6) 



Eq. (Q) generalizes Eq. What we have achieved 

here is that we manage to construct a large class of spin 
wave functions. Those spin wave functions $ are related, 
through the projection, to the mean-field wave function 
^ which have a simple form of Slater determinant. The 
constructed spin wave function can be labeled by the 
mean- field ansatz (uij,a l Q T l ). This allows us to study 
the quantum order of a spin wave function by studying 
the property of a simpler object, the mean-field ansatz 
(uij,a l T l ). 

In the study of phases and the related internal orders, 
the central question is to identify the universal properties 
of states. By definition, a universal property is a property 
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shared by all the states in the same phase. For classical 
systems, the symmetry of a state is a universal property 
We cannot change one state to another state without a 
phase transition if the two states have different symme- 
tries. Therefore, we can use the symmetry group (SG) 
to characterize the internal orders of classical states. 

Similarly, to characterize quantum orders, we need to 
find the universal properties of spin liquid wave function. 
Using the above projective construction, we can simplify 
the problem by considering the universal properties of the 
mean-field ansatz (uij,a l T l ) instead. Motivated by the 
classical systems, here we would like to propose that the 
symmetry of the mean- field ansatz (uij,a l T l ) is a uni- 
versal property. The symmetry group of the ansatz will 
be called the projective symmetry group (PSG). Under 
our conjecture, the quantum orders in spin liquids can be 
characterized by PSG's. 

At the first sight, it appears that the PSG character- 
ization of quantum orders is identical to the SG charac- 
terization of classical order used in Landau's theory In 
fact the PSG of a spin wave function is different from the 
SG of that wave function and the two characterizations 
are different. This is because the (uij,a T ) labeling of 
the physical spin wave function ' a o r! ) ig not an one- 
to-one labeling. Two mean-field ansatz (uij,a l r l ) and 
(uij, a l T l ) differed by an SU(2) gauge transformation de- 
fined by 

Uij = G(i) Uij G\j), a' (i)T l = G(i)a l (i)T l GHi), 
fa = G{i)i>i, G(i) e 5(7(2) (7) 

give rise to the same spin wave function ^("•j^o 7 "') = 
Q(u i:j ,a T ) "phis is because the up-spin (the boson), as 
a bound state of and ip 2 , is an singlet of the above 
gauge SU{2). Thus the spin wave function, as a pro- 
jected mean-field wave function, is invariant under the 
local SU{2) gauge transformation Eq. (Q). 

Due to the many-to-one labeling, an interesting situa- 
tion appears. In order for a spin wave function to have 
a translation symmetry, its corresponding ansatz is only 
required to be translation invariant up to an SU (2) gauge 
transformation. That is the ansatz should be invariant 
under translation followed by a proper gauge transforma- 
tion. For example, to have a translation symmetry in x 
direction, should satisfy 

Uij = G x T x ( Uij ), a l (i)T l = G x T x (a l (i)T l ); 
T x {uij) = Ui^,j-x, G x (uij) ee G x (i)uijG x (j) (8) 

for a certain SU(2) gauge transformation G x (i). For two 
spin wave functions with the same translation symmetry, 
their ansatz can be invariant under translation followed 
by different gauge transformations. Thus two spin liq- 
uids with the same symmetry can have different PSG's. 
This indicates that the PSG characterization is more re- 
fined then the SG characterization. PSG can describe 
those internal structures that cannot be distinguished by 
SG. Therefore, we can use PSG to characterize quantum 



orders which cannot be completely characterized by sym- 
metries. 

Now let us explicitly write down PSG for some simple 
ansatz. We will consider spin liquids with only transla- 
tion symmetry. The SG is generated by two translations 
in x and y directions SG — {T x ,T y }. The first ansatz 
is called Z2A ansatz which has a form it^i+m = u mi 
where u m = u < _ rn are generic 2x2 complex matri- 
ces. An element in PSG, in general, is formed by the 
combined transformation as in Eq. (|8j). Including the 
translation in the y-dircction, we find the PSG of the 
Z2A ansatz is generated by the following transforma- 
tions {Go, G X T X , G y Ty}, where three gauge transforma- 
tions Go, G x and G y are given by Go(i) = — r°, G x (i) = 
t°, G y {i) = t°. Here r° is the 2x2 identity matrix. 
Since the ansatz is already invariant under T x and T y , 
hence G x and G y are trivial. 

We would like to point out that a PSG contains a spe- 
cial subgroup, which will be called the invariant gauge 
group (IGG). An IGG is formed by pure gauge transfor- 
mations that leave the ansatz unchanged 

IGG={G\uij = G{i)uijG\j), (9) 
a l (i)r l = G(i)a l (i)r l GHi)} 

One can show that PSG, IGG, and SG are related 
PSG/ IGG — SG, thus comes the name projective sym- 
metry group. For the Z2A ansatz, we find the IGG is a 
Z 2 group generated by Go = — r°. Because of this we 
will call such a spin liquid a Z 2 spin liquid. 

Next we consider another ansatz, the Z2B ansatz, 
UiA +m = {-l) m ^u m . The Z2B PSG is still generated 
by {G ,G x T x ,GyT y }, but with a different G y : G {i) = 
-t°, G x (i) = t°, G y (i) = (-1) 1 -t°. Under transla- 
tion Ty, U iti + m — > Ui-y ti + m -y = ( ~ ) "** U^i+m • TllUS WC 

need a nontrivial gauge transformation G y — (— ) l * to 
remove the extra factor (— )™ la! . The IGG for the Z2B 
ansatz, generated by Go, is also Z 2 and the ansatz de- 
scribes another Z 2 spin liquid. 

One can show that[§ the above Z2A and Z2B PSG's 
cannot be transformed into each other by the SU(2) 
gauge transformation Eq. (^). Therefore the Z2A and 
Z2B ansatz describe two spin liquids with the same sym- 
metry but different quantum orders. This demonstrates 
that the ground state wave functions of spin liquids con- 
tain structures that cannot be characterized by symme- 
try. A generalization of the SG, PSG, can capture some 
of those extra structures. 

To experimentally measure the different quantum or- 
ders in the Z2A and Z2B states, we can measure the 
spectrum of spin-1 excitations [ie the two-spinon excita- 
tions). One can show that (see Eq. (96) in Ref. p|), for 
the Z2B state, the spin-1 spectrum is periodic in fc-space 
with a period 7r in both x and y directions. For the Z2A 
state, the spin-1 spectrum has the usual period of 2tt. 

After constructing the spin wave functions of some 
spin liquids via the projection Eq. @, we would like 
to ask: do those spin liquids actually exist for spin-1/2 
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systems? Are there any Hamiltonians such that the con- 
structed spin wave functions are the ground state of those 
Hamiltonians? This is a very hard question and there 
are some recent work supporting the existence of spin 
liquids. ]l3|, pH |l5| , It is also known that spin liquids 
do exist for certain large N systems. [0 |8|,]rJ. Here 
we would like to address an easier question about the 
stability of spin liquids: can a small perturbation in a 
spin Hamiltonian change the quantum order in the cor- 
responding spin liquid. This question can be addressed 
without knowing the details of the lattice spin Hamilto- 
nian. Wc only need to know the low lying spin excitations 
and their interactions. 

As an example, we consider the stability of the Z2A 
spin liquid. Within the SU(2) slave-boson mean-field 
theory, || |j| the spin excitations are described by the free 
spin- 1/2 spinons tp in the mean-field Hamiltonian Eq. (|^). 
In general, if Tr(u m )| is much less then \Tr(u m T l )\, the 
spinons have finite energy gap. As we go beyond the 
mean-field approximation, the low energy excitation also 
contain a collective mode described by the fluctuations 
Siiij. The spinons are no longer free since they interact 
with Suij. As pointed out in Ref. ^0|, the fluctuations 
Suij correspond to a gauge field at low energies. The 
gauge group of the low energy gauge fluctuations is deter- 
mined by the PSG.J01 In fact it is the IGG of the ansatz. 
(See also Ref. ||, Thus the Z 2 spin liquid has 

a Z 2 gauge fluctuations at low energies. The Z 2 gauge 
fluctuations will cause an interaction between the spinons 
since they carry unit Z 2 charge. However, the Z 2 gauge 
interaction between the spinons is short ranged. Thus 
even beyond the mean-field theory, the spin 1/2 spinon 
excitations are free at low energies and the low energy 
gauge fluctuations do not cause any instability. All ex- 
citations, spinons and the Z 2 gauge fluctuations, have a 
finite energy gap. In this way, we showed that the Z2A 
spin liquid can be a stable spin liquid which represents 
a quantum phase. The low energy properties and the 
quantum order of the spin liquid do not change as we 
perturb the spin Hamiltonian. 

In the above, we only showed that fluctuations around 
the Z 2 mean-field state do not cause any infrared di- 
vergence. However, even short distance fluctuations can 
cause instability if they are strong enough. Indeed the 
short distance fluctuations in our model are of order 



0(1). This means that the projection Eq. (^) causes a 
big change and it is unclear if the projected state and the 
original state share similar physical properties. To make 
our approach here into a controlled calculation, we need 
to generalize our model to some large-TV model. One 
way to do so is to generalize our mean-field Hamiltonian 
Eq. (§) to 

ij i 

where / = 1 , . . . , N and we have N copies of two- 
component fields ipi. After integrating out the fermions, 
the effective Lagrangian for the fluctuations Suij has a 
form C = NCo(Suij). In the large N limit, the fluctua- 
tions are suppressed and the mean-field theory becomes 
exact. Thus the two Z 2 states discussed above should 
exist in those large- N systems. 

To obtain the physical system that corresponds to 
the mean-field theory Eq. ([lO]), we note that a physi- 
cal state must be a singlet of the gauge SU(2) on every 
site. Thus the physical states on each site is formed by 
the singlet state of the fermions xpj. They include |0), 
ip* a jipp j£ Q ^|0), etc . The total number of physical state 

on each site is N p = Y^lo 2 JV ~ 3m (2m)!/m!. Such a 
system is a system with N p state per site (which can be 
viewed as a spin S = (N p — l)/2 system without any spin 
rotation symmetry). 

The PSG not only provides a concrete description of 
quantum order, it also allows us to partially classify quan- 
tum orders. || To understand the significance of PSG, we 
note that our understanding of solids is built on two cor- 
ner stones: (A) Solids contain an order that is related 
to broken symmetry. (B) The order can be described 
and classified by SG's. In this paper we propose that a 
understanding of quantum phase should be built on (at 
least) the following two concepts: (A) A quantum phase 
in general contain a new kind of order - quantum order 
- which may not have broken symmetry and local order 
parameters. (B) The quantum order can be (partially) 
described and classified by PSG's. 

This research is supported by NSF Grant No. DMR- 
97-14198. 
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